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Abstract. Weakly Interacting Massive Particle (WIMP) direct detec- 
tion experiments are just reaching the sensitivity required to detect Galac- 
tic dark matter in the form of neutralinos (or indeed any stable weakly 
interacting particle). Detection strategies and data analyses are often 
based on the simplifying assumption of a standard spherical, isothermal 
halo model, but observations and numerical simulations indicate that 
galaxy halos are in fact triaxial and anisotropic, and contain substruc- 
ture. The annual modulation and direction dependence of the event rate 
(due to the motion of the Earth) provide the best prospects of distin- 
guishing WIMP scattering from background events, however these signals 
depend sensitively on the local WIMP velocity distribution. I briefly re- 
view the status of WIMP direct detection experiments before discussing 
the dependence of the annual modulation signal on astrophysical input, 
in particular the structure of the Milky Way halo, and the possibility that 
the local WIMP distribution is not smooth. 



1. WIMPs a (very) brief introduction 

Any stable weakly interacting massive particle (WIMP) in thermal equilibrium 
in the early universe will generically have an interesting present day density, 
^wimp ~ O(^cdm) ~ 0.3. Furthermore supersymmetry provides a natural 
WIMP candidate, the lightest supersymmetric particle, the neutralino. There 
are basically two methods of detecting WIMPs: indirect detection, which in- 
volves detecting the products of WIMP annihilation (7, v, p, e + ), and direct 
detection, which involves detecting the energy deposited in a detector due to 
elastic scattering of WIMPs on the detector nuclei. I will focus on WIMP direct 
detection. For a review of particle dark matter see e.g. Bergstrom (2000). 

2. Direct detection signals 

Direct detection experiments are just reaching the sensitivity required to de- 
tect WIMPs. The expected event rates are very small ( O(10" 5 - 10) counts 
kg _1 day _1 ) and distinguishing a putative Weakly Interacting Massive Particle 
(WIMP) signal from backgrounds, such as neutrons from cosmic-ray induced 
muons or natural radioactivity, is crucial. The event rate depends on the veloc- 
ity of the detector relative to the Galactic rest frame and the Earth's motion (as 
shown in Fig. 1) provides two potential WIMP smoking guns. Firstly the event 
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Figure 1. A sketch of the motion of the Earth. The Sun moves about 
the Galactic centre with circular velocity v c ~ 220 km s -1 , while the 
Earth orbits the Sun with speed v e = 30km s" 1 in a close to circular 

o 

orbit, inclined at an angle of roughly 60 to the Galactic plane. 



rate is direction dependent, being greatest in the forward direction (Spergel 
1988). Secondly the Earth's velocity, and hence the event rate, varies annually 
due to the Earth's orbit about the Sun (Drukier, Freese & Spergel 1986). If 
the local WIMP velocity distribution is isotropic then the annual modulation 
is roughly sinusoidal with a maximum in early June (when the Earth's speed 
with respect to the Galactic rest frame is largest) and amplitude of order a few 
per-cent. 

There are currently more than 20 WIMP direct detection experiments be- 
ing carried out around the world. I will focus on those currently producing the 
most interesting results. The DAMA collaboration, using Nal with an exposure 
of 108 000 kg-days at Gran Sasso, have detected an annual modulation with 
the properties described above which they interpret as a WIMP signal (Bern- 
abei et al. 2003). Assuming a standard halo model with a Maxwellian velocity 
distribution with dispersion a = 270km s _1 (corresponding to an asymptotic 
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halo circular velocity of v c = 220 km s -1 ), they find a best fit WIMP mass 
m x £3 50 GeV and scattering cross-section (it«7x 10 -46 m 2 1 . Also of inter- 
est are three other experiments using different targets and a different strategy: 
Cryogenic Dark Matter Search (Ge target, 28 kd-days exposure at the Stan- 
ford Underground Facility (Akerib et al. 2003)), Edelweiss (Ge, 32 kg-days at 
Modane (Benoit et al. 2002)) and Zeplin I (liquid Xe, 230 kg-days at Boulby 
(Barton et al. 2002)). With their smaller exposures these experiments aim to 
constrain the mean WIMP scattering rate, rather than attempting to detect the 
annual modulation signal. The full DAMA allowed region of WIMP mass cross- 
section parameter space is shown in fig. 2, along with the exclusion limits from 
the CDMS, Edelweiss and Zeplin I experiments. The allowed region and the ex- 
clusion limits are all calculated assuming the standard halo model as described 
above and, taken at face value, the allowed region appears to be incompatible 
with the exclusion limits. 



3. Local velocity distribution 

The differential elastic scattering rate depends on the local WIMP density, p x , 
and the normalised WIMP speed distribution, in the rest frame of the detector, 



where v m { n is the minimum WIMP velocity which can kinematically produce a 
nuclear recoil of energy E. 

As discussed above, data analyses usually assume that WIMPs have an 
isotropic Maxwellian velocity distribution (i.e. that the Milky Way halo is an 
isothermal sphere), but observed and simulated halos are triaxial, anisotropic 
and contain substructure. Exclusion limits, which depend on the time aver- 
aged speed distribution, vary by of order tens of per-cent when triaxial and 
anisotropic halo models (with parameters chosen to match the properties of 
observed and simulated halos) are considered, with the shift being experiment 
dependent (Green 2002). 

The annual modulation signal is far more sensitive to the WIMP velocity 
distribution, and Belli et al. (2002) found that considering "non-standard" halo 
models led to a large increase in the size of the region of m x — a parameter space 
consistent with the DAMA annual modulation signal. Accurate calculation of 
the shape and phase of the annual modulation signal requires all three compo- 
nents of the Earth's velocity with respect to the Sun, and also the Sun's motion 
with respect to the Local Standard of Rest, to be taken into account (Green 
2003). Significantly, if the velocity distribution is not close to isotropic then 
the phase and shape of the annual modulation change and become incompatible 
with the DAMA annual modulation signal (Copi & Krauss 2003, Fornengo Sz 
Scopel 2003, Green 2003). 



x £ = p x /(0.3 GeV cm 3 ) parameterizes the uncertainty in the local WIMP density p x . 
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DATA listed top to bottom on plot 
CDMS June 2003, bkgd subtracted 

DAMA 2000 58k kg-days Nal Ann.Mod. 3sigma,w/o DAMA 1996 limit 

ZEPLIN I Preliminary 2002 result 

Edelweiss, 32 kg-days Ge 2000+2002+2003 limit 

Edelweiss 2 projection 

CDMS, projected at Soudan mine 

ZEPLIN 4 projection 
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Figure 2. Selected experimental results, plotted using the interactive 
limit plotter at http://dmtools.berkeley.edu/limitplots, which assumes 
the standard halo model. The solid region is the allowed region corre- 
sponding to the DAMA annual modulation signal, the solid lines the 
current exclusion limits from the CDMS, Edelweiss and Zeplin I exper- 
iments and the dotted lines the projected future sensitivities of these 
experiments. 
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4. (very) Small scale structure 

The velocity distributions of analytic halo models are derived via the Jean's 
equations, which assume that the phase space distribution function has reached 
a steady state. In Cold Dark Matter (CDM) cosmologies structure forms hier- 
archically, with galaxy halos forming from the merger and accretion of smaller 
subhalos (which themselves formed from even smaller subhalos), and the local 
velocity distribution may not have reached a steady state. Helmi, White &i 
Springel (2001) examined the velocity distribution of particles within a 4 kpc 
box located 8kpc from the centre of a simulated Milky Way like halo and found 
that, apart from a stream of fast moving particles from a late accreting subhalo, 
the velocity distribution was well approximated by a multi-variate gaussian. The 
WIMP direct detection event rate, however, depends on the dark matter dis- 
tribution on sub-mpc scales, many orders of magnitude smaller than the scales 
probed by even the highest resolution galaxy simulations. Moore et al. (2001) 
argued that the local velocity distribution will depend sensitively on the struc- 
ture and merger history of the first halos to form, while Stiff & Widrow (2003) 
have used a novel reverse technique to probe the velocity distribution at a single 
point in a simulation, finding that it appears to consist of a series of discrete 
peaks. 

The power spectrum on very small scales is an essential input for any at- 
tempt to study the fate of the first subhalos to form and the dark matter dis- 
tribution on very small scales. For neutralino CDM collisional damping and 
free-streaming erase power on comoving scales smaller than of order a pc (Hof- 
mann, Schwarz & Stocker 2001, Green, Hofmann & Schwarz 2003), and the first 
subhalos to form have mass of order 1O" 6 M (Green et al. 2003), 12 orders of 
magnitude lighter than the smallest subhalos resolved in galaxy simulations. 

Streams of WIMPs with small velocity dispersion (from late accreting sub- 
halos which pass through the solar neighbourhood or from the remnants of the 
first dense subhalos to form, should they survive to the present day) will pro- 
duce steps in the differential event rate, the position and amplitude of which 
vary annually. 



5. Summary 

Direct detection of WIMPs would confirm the existence of Cold Dark Matter 
(and probe particle physics beyond the standard model). Accurate astrophysical 
input (not just the local WIMP velocity distribution, but also the motion of the 
detector with respect to the Galactic rest frame) is required when calculating 
the WIMP annual modulation signal. Analyzing data assuming a sinusoidal 
modulation with fixed phase could lead to erroneous constraints on, or best fit 
values, for the WIMP mass and cross-section, even worse a WIMP signal could 
be overlooked. On the other hand using unrealistic halo models or parameter 
values could lead to overly restrictive exclusion limits or a misleadingly large 
range of allowed values of the WIMP mass and cross-section. Finally if WIMPs 
are directly detected then we will be able to probe the local velocity distribution 
and perhaps learn about the (sub-)structure of the Milky Way halo. 
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